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Spin-spin correlation functions of spin systems coupled to 2-d quantum 
gravity for < c < 1. 
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We perform Monte Carlo simulations of 2-d dynamically triangulated surfaces coupled to Ising and three-states 
Potts model matter. By measuring spin-spin correlation functions as a function of the geodesic distance we provide 
substantial evidence for a diverging correlation length at j3 c . The corresponding scaling exponents are directly 
related to the KPZ exponents of the matter fields as conjectured in H. 



1. INTRODUCTION 

The calculation of the dressed scaling ex- 
ponents is a milestone in the theory of two- 
dimensional quantum gravity (l]]. Strictly speak- 
ing, however, the derivation uses only finite size 
scaling arguments and involves only matter field 
correlators integrated over all space-time. The 
concept of a diverging correlation length at the 
critical point associated with the correlators is 
not used in the derivation and was never shown 
to actually exist. It is possible to define a correla- 
tion length by defining the two point correlator of 
a field <f>(x) in terms of the reparametrization in- 
variant geodesic distance R between two marked 
points: 

G (i?; A) = J V\g]V[<t>] e~ SG - SM J d 2 xd 2 y 
xV9(x)g(y) (Kx)(t>(y)5(D g (x,y) - R). (1) 

Even with this definition at hand, it is possible 
that the correlation length of the system coupled 
to gravity does not diverge, even though it does 
in flat space and the transition when coupled to 
gravity is continuous Q . The fixed volume corre- 
lator Gj,(R; V) is obtained by Laplace transform- 
ing G ( f,(R; A) and from it one obtains the correla- 
tors S^(R;V) = G^(R]V)/V Z(V), where Z(V) 



is the fixed volume partition function. The case 
<fi = 1 corresponds to the volume volume correla- 
tors Gi(i?;U), Gi(i?;A) and the average "area" 
of a spherical shell Si(R; V) or radius R. In the 
case of a matter field (f> with scaling dimension Ao 
in flat space, J °° dR S^R; V) ~ V^ A where A 
is the dressed scaling dimension of the field. From 
this scaling and the corresponding behaviour of 
S ( °\R;V) in flat space 



sf(i?;U) 



R 



R 2A 

we may conjecture that 
R d h -i 



f(x) = V^-^F^ix), (2) 



S4R;V) 



Rd h 



f( x ) = V^-^F^x), (3) 



where /(0) > 0, F(x) ~ x d h {i-A)-i for small 
x = R/V x / dh and d% is the fractal dimension of 
space-time as defined for example by the finite 
size scaling relation Si(R; V) ~ V 1 ~ 1 / dh F\{x). In 
this article we report on extensive numerical sim- 
ulations which provide substantial evidence that 
Eq. (|||) holds for matter fields with < c < 1. 
This points to the existence of a diverging cor- 
relation length of the conformal field theory cou- 
pled to gravity if we use the geodesic distance as 
a measure of length. 
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(a) Ising model on a fixed lattice 




Fit: c=3/4(fixed), C=4.224(4), B=-0.0044(l) 
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Figure 1. Data for F^ lat (x) as denned in Eq.(§) 
for small values of x = r/N^ dh , d h = 2. The fit 
is to Eq.@. 



2. SIMULATIONS AND RESULTS 

We perform numerical simulations on dynami- 
cally triangulated surfaces with S 2 topology with 
tadpoles and self-energy diagrams. We place Ising 
and three-states Potts model spins on the ver- 
tices of the surface which interact with the spins 
of the neighbouring vertices. The Monte Carlo 
updating of the triangulations is performed by 
the so-called flip algorithm and the spins are up- 
dated by standard cluster algorithms. The flips 
are organised in "sweeps" which consist of ap- 
proximately Nl accepted flips where Nl is the 
number of links of the triangulated surface. Af- 
ter a sweep we update the spin system ||. The 
results presented in this paper cover system sizes 
from 16000 to 128000 triangles and the number 
of sweeps is 1.7-5.0 x 10 6 . The simulations are 
performed at the analytically known infinite vol- 
ume (3 C Q. Geodesic distances r on the triangu- 
lations are defined as the shortest link distance 
between two vertices. The discretized volume of 
the system is the number of triangles Nt and we 
use the scaled distance x — R/N 1 ^ 11 , where N 
is the number of vertices. We measure the dis- 
cretized distributions corresponding to Si (R; V) 
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(a) Ising model coupled to gravity 






Fit: c=5/3 (fixed), C=0.945(4), B=0.0055{2) 
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(b) 3-state Potts model coupled to gravity 
Fit: c=7/5(fixed), C=1.045(l) 








o 


16000 




□ 


32000 




<> 


64000 




A 


128000 



0.05 0.1 0.15 0.2 0.25 0.3 



Figure 2. (a) Same as in Fig. for F^ sin9 (x) for 
4 =4. (6) Same as in (a) for F[ otts (x). 



and Sa,(R;V) 

m(r;N) = (J^CAi-r)), (4a) 

3 

n^r; N) = a t a 3 5(D i} - r)) , (4b) 

3 

where the indices i and j label vertices: i is a 
random fixed vertex for each measurement and j 
runs over all vertices of the given configuration. 
Dij is the link distance between the vertices la- 
belled by i and j. We determine the fractal di- 
mension dh by using ni(r;N) — N 1 ~ 1 / dh F{x) 
and n^r-N) = N^-^F^x). Our data is 
in very good agreement with the results obtained 
in ||,(|], namely dh = 4 for both the Ising model 
and the three-states Potts model coupled to grav- 
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Table 1 

The parameters of the fits to Eq. (^) for the Ising 
model on a flat T 2 for N T = 64000 (A = 1/8), 
and for the Ising model (A = 1/3) and the three- 
states Potts model (A = 2/5) on S 2 coupled to 
quantum gravity for Nt = 128000. 





A = 1/8 


A = 1/3 


A = 2/5 


Ci 


4.224(4) 


0.945(4) 


1.044(4) 


Bi 


-0.0044(1) 


0.0054(2) 


0.0003(3) 


c 


0.765(1) 


1.535(6) 


1.396(5) 


C 2 


4.39(2) 


0.834(7) 


1.040(8) 


c 


0.753(4) 


1.596(6) 


1.470(8) 


c 3 


4.25(4) 


0.870(7) 


1.10(1) 


B 3 


-0.0038(7) 


0.0036(1) 


0.0055(3) 



ity. 

We test the hypothesis of Eq. (||) by using 
the values of the dressed scaling dimensions of 
the Ising and three-states Potts model spin fields 
A = 1/3, 2/5. The prediction for d h = 4 is 



Ising 



Pottt 



(x) oc x 5/3 for (5a) 
(a;) oc x 7/5 for x <C 1 . (5b) 



In order to calibrate the expected accuracy with 
which one can determine the exponent d/j(l — 
A) — 1 we have performed simulations of an Ising 
model on a fixed triangular lattice with peri- 
odic boundary conditions (A = 1/8). We mea- 
sured F^ lat (x) of Eq. (|) and we expect that 
Fp at {x) oc a; 3 / 4 for x < 1. 

In Table || we show the results of the fits for the 
largest lattice to the following functional forms: 



Ft(x) = C lX d ^-^- 1 
= C 2 x c , 
= C 3 x c + B 3 . 



Bi, 



(6a) 
(6b) 
(6c) 



lating the largest lattices, 128000 triangles, is im- 
portant for obtaining enough data points in the 
relevant region x < 0.45. The fits to the pre- 
dicted behaviour Eq. ( |6a| ) are good and the finite 
size correction B approaches zero as the volume 
increases. Even if the exponent c is allowed to 
vary, as in the fits to Eqs. (6b) and (|6c|), it ap- 
proaches its predicted value convincingly as the 
volume is increased. The small discrepancy is 
consistent with finite size effects, which clearly are 
more important than the statistical errors quoted 
in Table ||. We find that the difference in the val- 
ues of c obtained from the fits to Eqs. (6b) and 
( |6c| ) gives a measure of the systematic errors en- 
tering from varying the range and type of the fits. 
For a more complete presentation of the data we 
refer the reader to . 
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